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The longitudinal and transversal spin decoherence times, T\ and T%, in semiconductor quantum 
dots are investigated from equation-of-motion approach for different magnetic fields, quantum dot 
sizes, and temperatures. Various mechanisms, such as the hyperfine interaction with the surrounding 
nuclei, the Dresselhaus spin-orbit coupling together with the electron-bulk-phonon interaction, the 
g-factor fluctuations, the direct spin-phonon coupling due to the phonon-induced strain, and the 
coaction of the electron-bulk/surface-phonon interaction together with the hyperfine interaction 
are included. The relative contributions from these spin decoherence mechanisms are compared in 
detail. In our calculation, the spin-orbit coupling is included in each mechanism and is shown to 
have marked effect in most cases. The equation-of-motion approach is applied in studying both the 
spin relaxation time T\ and the spin dephasing time T% % either in Markovian or in non-Markovian 
limit. When many levels are involved at finite temperature, we demonstrate how to obtain the 
spin relaxation time from the Fermi Golden rule in the limit of weak spin-orbit coupling. However, 
at high temperature and/or for large spin-orbit coupling, one has to use the equation-of-motion 
approach when many levels are involved. Moreover, spin dephasing can be much more efficient 
than spin relaxation at high temperature, though the two only differs by a factor of two at low 
temperature. 

PACS numbers: 72.25.Rb, 73.21.La,71.70.Ej 



I. INTRODUCTION 

One of the most important issues in the growing 
field of spintronics is quantum information processing 
in quantum dots (QDs) using electron spini^^ A 
main obstacle is that the electron spin is unavoidably 
coupled to the environment (such as, the lattice) which 
leads to considerable spin decoherence (including lon- 
gitudinal and transversal spin decoherences)^ 7 - Vari- 
ous mechanisms, such as, the hyperfine interaction with 
the surrounding nucleif^^ the Dresselhaus/Rashba spin- 
orbit coupling (SOC) 10 : 11 together with the electron- 
phonon interaction, <?-factor fluctuations^ the direct 
spin-phonon coupling due to the phonon-induced strain, 9 
and the coaction of the hyperfine interaction and the 
electron-phonon interaction can lead to the spin deco- 
herence. There are quite a lot of theoretical works 
on spin decoherence in QD. Specifically, Khaetskii and 
Nazarov analyzed the spin-flip transition rate using a 
perturbative approach due to the SOC together with the 
electron-phonon interaction, g- factor fluctuations, the di- 
rect spin-phonon coupling due to the phonon-induced 
strain qualitatively ! 13 : 14 : 15 After that, the longitudinal 
spin decoherence time T\ due to the Dresslhaus and/or 
the Rashba SOC together with the electron-phonon in- 
teraction were studied quantitatively in Refs . 1 1 6lJ 1 7lJ 1 8lJ 1 9l . 
HimiUHHIiiil Among these works, Cheng et al& 
developed an exact diagonalization method and showed 
that due to the strong SOC, the previous perturbation 
metho d 14 ! 15 ' 16 is inadequate in describing T±. Further- 
more, they also showed that, the perturbation method 



previously used missed an important second-order en- 
ergy correction and would yield qualitatively wrong re- 
sults if the energy correction is correctly included and 
only the lowest few states are kept as those in Refs. 

These results were later confirmed by Deste- 
fani and Ulloa^I The contribution of the coaction of the 
hyperfine interaction and the electron-phonon interac- 
tion to longitudinal spin decoherence was calculated in 
Refs. and [H. In contrast to the longitudinal spin de- 
coherence time, there are relatively fewer works on the 
transversal spin decoherence time, T2, also referred to as 
the spin dephasing time (while the longitudinal spin de- 
coherence time is referred to as the spin relaxation time 
for short). The spin dephasing time due to the Dressel- 
haus and/or the Rashba SOC together with the electron- 
phonon interaction was studied by Scmcnov and Kim 29 
and by Golovach et The contributions of the hyper- 
fine interaction and the q-factor fluctuation were studied 
in Refs J^pllllll^ and 
in Ref. [Ifl respectively. However, a quantitative calcula- 
tion of electron spin decoherence induced by the direct 
spin-phonon coupling due to phonon-induced strain in 
QDs is still missing. This is one of the issues we are 
going to present in this paper. In brief, the spin re- 
laxation/dephasing due to various mechanisms has been 
studied previously in many theoretical works. However, 
almost all of these works only focus individually on one 
mechanism. Khaetskii and Nazarov discussed the ef- 
fects of different mechanisms on the spin relaxation time. 
Nevertheless, their results are only qualitative and there 
is no comparison of the relative importance of the dif- 
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ferent mechanisms ! 13 ' 1 i 15 Recently, Semenov and Kim 
discussed various mechanisms contributed to the spin 
dephasing— where they gave a "phase diagram" to in- 
dicate the most important spin dephasing mechanism in 
Si QD where the SOC is not important. However, the 
SOC is very important in GaAs QDs. To fully under- 
stand the microscopic mechanisms of spin relaxation and 
dephasing, and to achieve control over the spin coherence 
in QDs ; 47 ' 48 ' 49 one needs to gain insight into the relative 
importance of each mechanism to T\ and T 2 under vari- 
ous conditions. This is one of the main purposes of this 
paper. 

Another issue we are going to address relates to differ- 
ent approaches used in the study of the spin relaxation 
time. The Fermi- Golden-rule approach, which is widely 
used in the literature, can be used in calculation of the re- 
laxation time Ti^f between any initial state \i) and final 

State [J) . 12.13.14.15.16.17.18.19.21.23.24.25.27.28.50.51.52 JJ QW _ 

ever, the problem is that when the process of the spin 
relaxation relates to many states, (e.g., when tempera- 
ture is high, the electron can distribute over many states), 
one should find a proper way to average over the relax- 
ation times (r^f) of the involved processes to give the 
total spin relaxation time (Ti). What makes it difficult 
in GaAs QDs, is that all the states are impure spin states 
with different expectation values of spin. In the existing 
literature, spin relaxation time is given by the average of 
the relaxation times of processes from the initial state \i) 
to the final state |/) (with opposite majority spin of \i)) 
weighted by the distribution of the initial states /i , 18 ' 51 ' 52 
i.e., 

IT* = C 1 ) 

if 

This is a good approximation in the limit of small SOC 
as each state only carries a small amount of minority 
spin. However, when the SOC is very strong which hap- 
pens at high levels, it is difficult to find the proper way 
to perform the average. We will show that Eq. (fl]) is 
not adequate any more. Thus, to investigate both T\ 
and T2 at finite temperature for arbitrary strength of 
SOC, we develop an equation-of-motion approach for the 
many-level system via projection operator technique^ 6 - in 
the Born approximation. With the rotating wave ap- 
proximation, we obtain a formal solution to the equation 
of motion. By assuming a proper initial distribution, 
we can calculate the evolution of the expectation value 
of spin. We thus obtain the spin relaxation/dephasing 
time by the 1/e decay of the expectation value of spin 
operator (S z ) or (to its equilibrium value), with 

S+ = S x + iS y . With this approach, we are able to study 
spin relaxation/dephasing for various temperature, SOC 
strength, and magnetic field. 

For quantum information processing based on electron 
spin in QDs, the quantum phase coherence is very im- 
portant. Thus, the spin dephasing time is a more rel- 
evant quantity. There are two kinds of spin dephasing 
times: the ensemble spin dephasing time TjT an d the 



irreversible spin dephasing time T2. For a direct mea- 
surement of an ensemble of QDs^ or an average over 
many measurements at different times where the config- 
urations of the environment have been changed ) 59 ' 60 ' 61 
it gives the ensemble spin dephasing time T 2 *. The ir- 
reversible spin dephasing time T2 can be obtained by 
spin echo measurement ! 60 ' 61 A widely discussed source 
which leads to both T 2 * an d T% is the hyperfine interac- 
tion between the electron spin and the nuclear spins of 
the lattice. As it has been found that T| is around 10 ns, 
which is too short and makes a practical quantum infor- 
mation processing difficult in electron spin based qubits 
in QDs. Thus a spin echo technique is needed to remove 
the free induction decay and to elongate the spin dephas- 
ing time. Fortunately, this technique has been achieved 
first by Petta et al. for two electron triplet-singlet system 
and then by Koppens et al. for a single electron spin sys- 
tem. The achieved spin dephasing time is ~ 1 /is, which 
is much longer than T 2 * . We therefore discuss only the ir- 
reversible spin dephasing time T2 throughout the paper, 
i.e., we do not consider the free induction decay in the 
hyperfine-interaction-induced spin dephasing. 

It is further noticed that Golovach et al. have shown 
that the spin dephasing time T 2 is two times the spin 
relaxation time T\*& However, as temperature increases, 
this relation does not hold. Semenov and Kim on the 
other hand reported that the spin dephasing time is much 
smaller than the spin relaxation timei^ 9 - In this paper, 
we calculate the temperature dependence of the ratio of 
the spin relaxation time to the spin dephasing time and 
analyze the underlying physics. 

This paper is organized as follows: In Sec. II, we 
present our model and formalism of the equation-of- 
motion approach. We also briefly introduce all the spin 
decoherence mechanisms considered in our calculations. 
In Sec. Ill we present our numerical results to indicate 
the contribution of each spin decoherence mechanism to 
spin relaxation/dephasing time under various conditions 
based on the equation-of-motion approach. Then we 
study the problem of how to obtain the spin relaxation 
time from the Fermi Golden rule when many levels arc 
involved in Sec. IV. The temperature dependence of T\ 
and T 2 is investigated in Sec. V. We conclude in Sec. VI. 



II. MODEL AND FORMALISM 

A. Model and Hamiltonian 

We consider a QD system, where the QD is confined 
by a parabolic potential V c (x, y) = ^m'w^i 2 + y 2 ) in 
the quantum well plane. The width of the quantum well 
is a. The external magnetic field B is along z direction, 
except in Sec. IV. The total Hamiltonian of the system 
of electron together with the lattice is: 

H T = H e + H L + H eL , (2) 
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where H e , Hl, H e L are the Hamiltonians of the elec- 
tron, the lattice and their interaction, respectively. The 
electron Hamiltonian is given by 

H e = ^ + V c {v)+H z +H S o (3) 

where P = -iW + f A with A = (B ± /2)(-y, x) (B±_ is 
the magnetic field along z direction), Hz = ^g^B^l ■ cr 
is the Zeeman energy with ps the Bohr magneton, 
and Hso is the Hamiltonian of SOC. In GaAs, when 
the quantum well width is small or the gate-voltage 
along the growth direction is small, the Rashba SOC is 
unimportant ^ Therefore, only the Dresselhaus term 1 ^ 
contributes to Hso- When the quantum well width is 
smaller than the QD radius, the dominant term in the 
Dresselhaus SOC reads 

Hso = ^(P^Xoi-P*^ + Py°y) , (4) 

with 70 denoting the Dresselhaus coefficient, Ao being 
the quantum well subband index of the lowest one and 
(P 2 )x = -h 2 J ip* zX (z)d 2 /dz 2 ip z x(z)dz. The Hamiltonian 
of the lattice consists of two parts Hl = H p h + H nuc i e i , 
where H p h = J2qr/ ^q^qr^qij (aV a i s the phonon cre- 
ation/annihilation operator) describes the vibration of 
the lattice and H nuc i ei = J2j 7/ B ' Ij (7/ is tne gyro- 
magnetic ratios of the nuclei and lj is the spin of the j- 
th nucleus) describes the precession of the nuclear spins 
of the lattice in the external magnetic field. We focus 
on the spin dynamics due to hyperfine interaction at a 
time scale much smaller than the nuclear dipole-dipole 
correlation time (10 -4 s in GaA a 33 i 40 ), where the nuclear 
dipole-dipole interaction can be ignored. Under this ap- 
proximation, the equation of motion for the reduced elec- 
tron system can be obtained which only depends on the 
initial distribution of the nuclear spin bathi^ 3 . The in- 
teraction between the electron and the lattice also has 
two parts H e L = H e i + H e - p h, where H e j is the hy- 
perfine interaction between the electron and nuclei and 
H e - p h represents the electron-phonon interaction which 
is further composed of the electron-bulk-phonon (BP) 
interaction H ep , the direct spin-phonon coupling due to 
the phonon-induced strain H str ain and phonon-induced 
g-factor fluctuation H g . 

B. Equation-of- motion approach 

The equations of motion can describe both the co- 
herent and the dissipative dynamics of the electron sys- 
tem. When the quasi-particles of the bath relax much 
faster than the electron system, the Markovian approx- 
imation can be made; otherwise the kinetics is the non- 
Markovian. For electron-phonon coupling, due to the 
fast relaxation of the phonon bath and the weak electron- 
phonon scattering, the kinetics of the electron is Marko- 
vian. Nevertheless, as the nuclear spin bath relaxes much 
slower than the electron spin, the kinetics due to the 



coupling with nuclei is of non- Markovian type^ -^^ It 
is further noted that there is also a contribution from 
the coaction of the electron-phonon and electron-nuclei 
couplings, which is a fourth order coupling to the bath. 
For this contribution, the decoherence of spin is mainly 
controlled by the electron-phonon scattering while the 
hyperfine (Overhauser) field 5 ^ acts as a static magnetic 
field. Thus, this fourth order coupling is also Markovian. 
Finally, since the electron orbit relaxation is much faster 
than the electron spin relaxation^ 5 - 5 - we always assume a 
thcrmo-cquilibrium initial distribution of the orbital de- 
grees of freedom. 

Generally, the interaction between the electron and the 
quasi-particle of the bath is weak. Therefore the first 
Born approximation is adequate in the treatment of the 
interaction. Under this approximation, the equation of 
motion for the electron system coupled to the lattice envi- 
ronment can be obtained with the help of the projection 
operator technique. 56 We then assume a sudden approxi- 
mation so that the initial distribution of the whole system 
is p(t = 0) = p e (0) ® p L {0), where p e , p L is the density 
matrix of the system and bath respectively. This approx- 
imation corresponds to a sudden injection of the electron 
into the quantum dot, which is reasonable for genuine 
experimental setupj^ 3 . As the initial distribution of the 
the lattice p L (0) commutates with the Hamiltonian of 
the lattice Hl , the equation of motion can be written as 

^ = - l -[H e + Tr L (H eL p L (0) ) , p e (t)] 

i r* 

-jp dTTi L [H eL ,U (T)(V [H eL ,p e {t-r) 

®p L m)Ul{r)] , (5) 

where p e {t) is the density operator of the electron sys- 
tem at time i, Tr^ stands for the trace over the lattice 
degree of freedom, and Uq(t) — q-^Hl+h^t - IS ^j me _ 

evolution operator without H e ^. V = 1 — p L {0) £g) Tr^ 
is the projection operator. The initial distribution of the 
phonon system is chosen to be the thermo-equilibrium 
distribution. 2 - It has been shown by previous theoretical 
studies that the initial state of the nuclear spin bath is 
crucial to the spin dephasing and relaxationi 30 i 32 i 33 Al- 
though it may take a long time (e.g., seconds) for the 
nuclear spin system to relax to its thermo-equilibrium 
state, one can still assume that its initial state is the 
thermo-equilibrium one. This assumption corresponds 
to the genuine case of enough long waiting time during 
every individual measurement. For a typical setup at 
above 10 mK and with about 10 T external magnetic 
field, the thermo-equilibrium distribution is a distribu- 
tion with equal probability on every state. For these ini- 
tial distributions of phonons and nuclear spins, the term 
Tr L (iJ eL p L (0)) is zero. Thus, 

V[H eL ,p e (t- T )®p L {0)] = [H eLlP e {t-T)®p L {0)} . (6) 
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The equation of motion is then simplified to, 



d P e {t) 
dt 



= -~[H e , p e (t)} - i J drTr L [H eL , [fl^(-r), 



US(t)p Ie (t-r)USHt)p L (0)}] , 



(7) 



where -ffg L and p Ie are the corresponding operators (B & l 
and // ) in the interaction picture, and U§(t) = e~ lH,!t 
is the time-evolution operator of H e . It should be fur- 
ther noted that the first Born approximation can not 
fully account for the non-Markovian dynamics due to the 
hyperfine interaction with nuclear spins i 33 ' 57 Only when 
the Zeeman splitting is much larger than the fluctuating 
Overhauser shift, the first Born approximation is ade- 
quate. For GaAs QDs, this requires B > 3.5 TM In this 
paper, we focus on the study of spin dephasing for the 
high magnetic field regime of B > 3.5 T under the first 
Born approximation, where the second Born approxima- 
tion only affects the long-time behavior.— Later we will 
argue that this correction of long time dynamics changes 
the spin dephasing time very little. 



1. Markovian kinetics 

The kinetics due to the coupling with phonons can be 
investigated within the Markovian approximation, where 
the equation of motion reduces to, 



i 1 /"' 

= -^[H e ,p e (t)}- dTTr ph [H e - ph , 

[H I e - ph (-T),p B (t)® f ? h (p)]]. (8) 



d P e (t) 
dt 



Here Tr p ^ is the trace over phonon degrees of freedom 
and p ph (0) is the initial distribution of the phonon bath. 
Within the basis of the eigen-states of the electron Hamil- 
tonian, {|^)}, the above equation reads, 



d e _ ; (gli ~ £ ll) e 



Put, ® pI 



-H, 



e 3 e 4 

I e—ph e 



Pli^P P e q H!J 2 h ) + H.c} . (9) 



Here H^f = (h\H^ ph \£ 3 ) and H^ ph = 

(£i\H* h (— t)\£s). A general form of the electron- 
phonon interaction reads 



H, 



$ qJ) (a q , + at )X qJ7 (r,o-) 



(10) 



Here, r\ represents the phonon branch index; is the 
matrix element of the electron-phonon interaction; a q?; 
is the phonon annihilation operator; A qj; (r, cr) denotes 
a function of electron position and spin. Substituting 



this into Eq. ([9]), we obtain, after integration within the 
Markovian approximation^ 



d e 



(6^ - ei 2 ) 



qq* e 



xC q ,( £f4 - £f3 )-^lS>^ 4 
xC q „(e, 3 -e tl )}+ B.c. } (11) 

in which X?X_ = (4|A q „(r,<r)|f 2 ), and C w (Ae) = 
n((J m ) S( Ae + w q „) + [n {w m ) + 1] 8( Ae - u> qr) ) . Here n (u m ) 
represents the Bose distribution function. Equation (jlip 
can be written in a more compact form 



i 3 ii 



(12) 



which is a linear differential equation. This equation can 
be solved by diagonalizing A. Given an initial distribu- 
tion p! ^ (0), the density matrix pf £ (t) and the expec- 
tation value of any physical quantity (0)t = Tr(Op e (i)) 
at time t can be obtained^ 

(0) t = Tr(Op e ) 

ti-te 

xrr,W4, ' p (4)(w/w 6 (°) ( 13 ) 

with r = P 1 AP being the diagonal matrix and P repre- 
senting the transformation matrix. To study spin dynam- 
ics, we calculate (S z ) t (\(S + ) t \) and define the spin relax- 
ation (dephasing) time as the time when (S z ) t (\(S + ) t \) 
decays to 1/e of its initial value (to its equilibrium value). 



2. Non-Markovian kinetics 

Experiments have already shown that for a large en- 
semble of quantum dots or for an ensemble of many mea- 
surements on the same quantum dot at different times, 
the spin dephasing time due to hyperfine interaction is 
quite short, ~ 10 ns i 58 ' 59 i 60 ' 61 This rapid spin dephas- 
ing is caused by the ensemble broadening of the preces- 
sion frequency due to the hyperfine fields4& When the 
external magnetic field is much larger than the random 
Overhauser field, the rotation due to the Overhauser field 
perpendicular to the magnetic field is blocked. Only 
the broadening of the Overhauser field parallel to the 
magnetic field contribute to the spin dephasing. To de- 
scribe this free induction decay for this high magnetic 
field case, we write the hyperfine interaction into two 
parts: B eI = h • S = B eI1 + B eI2 . Here h = (h x , h y , h z ) 
and S = (S x , S y , S z ) are the Overhauser field and the 
electron spin respectively. B e ji = h z S z and B e i2 = 
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7}(h + S- + h_S + ) with h± = h x ± ih y . The longitudi- 
nal part H e ji is responsible for the free induction de- 
cay, while the transversal part H e j2 is responsible for 
high order irreversible decay. As the rapid free induction 
decay can be removed by spin echo ; 60 ' 61 elongating the 
spin dephasing time to ~ 1 /is which is more favorable 
for quantum computation and quantum information pro- 
cessing, we then discuss only the irreversible decay. We 
first classify the states of the nuclear spin system with its 
polarization. Then wc reconstruct the states within the 
same class to make it spatially uniform. These uniformly 
polarized pure states, |n)'s, are eigen-states of h z . They 
also form a complete-orthogonal basis of the nuclear spin 
system. A formal expression of |n) iaM 



E 



N 



(14) 



Here \I,rrij} denotes the eigen-state of the z-component 
of the j-th nuclear spin Ij Z with the eigenvalue Hrrij . N 
denotes the number of the nuclei. The equation of motion 
for the case with initial nuclear spin state Pi s (0) = |n)(n| 
is given by^ 3 - 



dp e (t) 
dt 



1 /"* 

J dTTr ns [H eI2 ,Uf(T) 



x [H eI2 , p*(t - t) <g> pTmufir)] ■ (15) 

As in traditional projection operator technique, the dy- 
namics of the nuclear spin subsystem is incorporated 
self-consistently in the last ternii 33 i 56 Here Tr ns is the 
trace over nuclear spin degrees of freedom. Uq 1 (t) = 
exp[— ir{H e + Hj + H e n)]. The Ovcrhauser field is given 
by h = J2j AvoljS(r — Rj), where the constants A and 
vq are given later. Ij and R^ are the spin and posi- 
tion of j-th nucleus respectively. As mentioned above, 
the initial state of the nuclear spin bath is chosen to be 
a state with equal probability of each state, therefore 
P" s (0) = £ n l^»R with N w = E„l being the 
number of states of the basis To quantify the ir- 

reversible decay, we calculate the time evolution of 
for every case with initial nuclear spin state \n). We then 
sum over n and obtain 



ll<s+) t || = £l<£ 



(«)\ 



(16) 



It is noted that the summation is performed after the 
absolute value of (S^)t- Therefore, the destructive in- 
terference due to the difference in precession frequency 
ujzn-, which originates from the longitudinal part of the 
hyperfine interaction (H e n), is removed. We thus use 
1 /e decay of the envelope of 1 1 (S + ) t \ | to describe the irre- 
versible spin dephasing time T%. Similar description has 
been used in the irreversible spin dephasing in semicon- 
ductor quantum wells22 and the irreversible inter-band 
optical dephasing in semiconductors! 6 ^ 64 



Expanding Eq. (fT"5| in the basis of one obtains, 

£3 

J ^ _ n n 



ni l 3 t 



-Hti-ltAAt ~ ^XLmJ + H.c) . (17) 

Here H%. nita = (rd x \H eI2 \n x Zz) and H^l. = 
(n£i\H^ I2 (—T)\ni£ 3 ). For simplicity, we neglect the 
terms concerning different orbital wavefunctions which 
are much smaller. For small spin mixing, assuming an 
equilibrium distribution in orbital degree of freedom, un- 
der rotating wave approximation, and trace over the or- 
bital degree of freedom, we finally arrive at 



dt 



knn' 



X [/l— Jfcn'n + [h-]knn'[h+]kn'n) 

X exp[fr(u/ fcn - uJkn>)}}(S+ )t-r ■ (18) 

Here uj zn = J2k fk{E zk /h + oj kn ) with E zk representing 
the electron Zeeman splitting of the fc-th orbital level. 
[hi\knn> = (n\(k\hi\k)\n') (i = ±,z). u kn = [h z ] knn + e nz 
with e nz denoting the nuclear Zeeman splitting which is 
very small and can be neglected. By solving the above 
equation, we obtain | (S+^)t\ for a given \n) . We then sum 
over n and determine the irreversible spin dephasing time 
T 2 as 1 je decay of the envelop of 1 1 (S+)t| | ■ By noting that 
only the polarization of nuclear spin state \n) determines 

the evolution of \(S^)t\, the summation over n is then 
reduced to summation over polarization which becomes a 
integration for large N. This integration can be handled 
numerically. 

In the limiting case of zero SOC and very low tem- 
perature, only the lowest two Zeeman sublevels are con- 
cerned. The equation for (S+)t with initial nuclear spin 
state Pi s (0) = |n) reduce to 

— {S + )t=iu} m {S + )t-^ J <M~£([/i+]„„/ 

n' 

x [h-] n i n + [/i_] nn /[/i + ] n / n )exp[ir(a; n - uv)]}(S+) t _ r 



= iw z (S+)t- / dTY,(T)(S + )t-r 
Jo 



(19) 



In this equation uj zn = {gpsB + [h z ] nn >)/h, [h^] nn ' = 
(n\(ij}i\h^\i^i)\n') (£ = ±,z and ipi is the orbital quantum 
number of the ground state), and w„ = [h z ] nn . Similar 
equation has been obtained by Coish and Loss,— and 
later by Deng and Hu 3 -^ at very low temperature such 
that only the lowest two Zeeman sublevels are considered. 
Coish and Loss also presented an efficient way to evaluate 
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£(t) in terms of their Laplace transformations, S(s) 
/ °° dTe~ ST Y>{T). They gave, 



AK 

n' 

+ [h-]nn'[h+\n'n)/{s ~ i8uJ nn >) , (20) 

with 8uj n n' — \{u n — w n /). With the help of this tech- 
nique, we are able to investigate the spin dephasing due 
to the hyperfine interaction. 



C. Spin decoherence mechanisms 

In this subsection we briefly summarize all the spin de- 
coherence mechanisms. It is noted that the SOC modifies 
all the mechanisms. This is because the SOC modifies the 
Zeeman splitting^ and the spin-resolved eigen-states of 
the electron Hamiltonian, it hence greatly changes the 
effect of the electron-BP scattering^ These two modifi- 
cations, especially the modification of the Zeeman split- 
ting, also change the effect of other mechanisms, such 
as, the direct spin-phonon coupling due to the phonon- 
induced strain, the g-factor fluctuation, the coaction of 
the electron-phonon interaction and the hyperfine inter- 
action. In the literature, except for the electron-BP scat- 
tering, the effects from the SOC are neglected except for 
the work by Woods et al*£ in which the spin relaxation 
time between the two Zeeman sub-levels of the lowest 
electronic state due to the phonon-induced strain is in- 
vestigated. However, the perturbation method they used 
does not include the important second-order energy cor- 
rection. In our investigation, the effects of the SOC are 
included in all the mechanisms and we will show that 
they lead to marked effects in most cases. 



static dielectric constant. The acoustic phonon spectra 
ujqqi = v s iq for the longitudinal mode and ujq P t = v st q 
for the transverse mode with v s i and v st representing 
the corresponding sound velocities. 

Besides the electron-BP scattering, electron also cou- 
ples to vibrations of the confining potential, i.e., the 
surface-phonons ^ 



2puj c , v V 



(a, 



)e ql) • V r K(r) , 



(22) 

in which e ql) is the polarization vector of a phonon mode 
with wave-vector q in branch 77. However, this contri- 
bution is much smaller than the electron-BP coupling. 
Compared to the coupling due to the deformation poten- 
tial for example, the ratio of the two coupling strengths is 
« huJofEqlo , where Iq is the characteristic length of the 
quantum dot and hu>o is the orbital level splitting. The 
phonon wave- vector q is determined by the energy differ- 
ence between the final and initial states of the transition. 
Typically phonon transitions between Zeeman sublevels 
and different orbital levels, qlo ranges from 0.1 to 10. 
Bearing in mind that hu>o is about 1 meV while 5 = 7 eV 
in GaAs, fkoo/Eql is about 10~ 3 . The piezoelectric cou- 
pling is of the same order as the deformation potential. 
Therefore spin decoherence due to the electron-surface- 
phonon coupling is negligible. 



2. Direct spin-phonon coupling due to phonon-induced 
strain 



The direct spin-phonon coupling due to the phonon- 
induced strain is given by^ 



1. SOC together with electron-phonon scattering 

As the SOC mixes different spins, the electron-BP scat- 
tering can induce spin relaxation and dephasing. The 
electron-BP coupling is given by 



Hep — ^ ^ -^qr; (^qr? 
q?f 



-qjj 



(21) 



where Af qr) is the matrix element of the electron-phonon 
interaction. In the general form of the electron phonon 
interaction H e - P hi $ qt? = M qt? and X q7) (r, <r) = e 4qr . 
l-ft^qszl 2 = hE 2 q/2pv s iV for the electron-BP coupling due 
to the deformation potential. For the piezoelectric cou- 
pling, |M qpi | 2 = (32fnT 2 e 2 e 2 14 /K, 2 pv s iV)[(3q x q y q z ) 2 /q 7 } 
for the longitudinal phonon mode and J2j=i 2 l-^qptj \ 2 = 
[i2h^e 2 e 2 J{n 2 pv st q 5 V)][q 2 x q 2 y + q 2 y q 2 z + q 2 z q\ 
(3q x q y q z ) 2 /q 2 ] for the two transverse modes. Here 
S stands for the acoustic deformation potential; p is the 
GaAs volume density; V is the volume of the lattice; 
ei4 is the piezoelectric constant and k denotes the 



H strain = ~h* '(p) ' C , 



where hi 



-Dp x (e yy - e zz ), hi 



-Dp y (e z 



(23) 
t xx ) and 

K = -Dp z (e xx -e yy ) with p = {p Xl p yi p z ) = -ihV and 
D being the material strain constant, ey — x,y,z) 
can be expressed by the phonon creation and annihilation 
operators: 



qn=l,t 1 ,t 2 



2 W 2puj cm V 



( a q,r; + a -q,n)(^irjqj 

+ ^K q ' r , (24) 



in which £n — qijq for the longitudinal phonon 
mode and (£ xtl , £ ytl , £ ztl ) = (q x q z ,q y q z ,-q 2 )/qq\\, 
(£,xt 2 ,£, y t 2 ,£zt 2 ) = {q y ,-qx,0)/q\\ for the two trans- 
verse phonon modes with qu = 



q 2 



q 2 . There- 



fore, in the general form of electron-phonon interaction 
H e -ph, $ qr , = -iD^/h/{Z2puj m V) and X m (r,a) = 
Y,i 3 k e ijk{£,]v1j - 6c»?<2fc)fte M1 ' r (7 l with e i]k denoting the 
Levi-Civita tensor. 
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3. g-factor fluctuation 

The spin-lattice interaction via phonon modulation of 
the (/-factor is given byi£ 



ijkl—x,y,z 



AijklHBBi<?j£kl ) 



(25) 



where e k i is given in Eq. {24J) and A^m is a 
tensor determined by the material. Therefore in 
H e - ph , $ qr , = iy / h/(32pujq, 1 V) and A q?) (r,er) = 
fc,i A id!k! HJ, B Bi(£ kri q k - £ it? q ; )cr. ) e lqr . Due to the 
axial symmetry with respect to the z-axis, and keep- 
ing in mind that the external magnetic field is along 
the z direction, the only finite element of H g is H g = 
{(A 33 -A 3 i)e zz +A 3 i J2i e zi ]hiJ,BBa z /2 with A 33 = A zzzz , 
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A z 



and A$q = A. 



*xyxy • 



A 33 + 2A 31 = 0. 



and 



VI 



(3) 



(£ 2 \H ep \m)(m\H eI (r)\£ 1 ) 



el—ph 



(£ 2 \H eI (r)\m){m\H ep \£ x ) 



E 



£i 2 — £ r. 



}(ti\ ,(29) 



in which \£i) and \£ 2 ) are the eigen states of H e . By using 
the notations in H e _ p h, $ q?? — j-y/h/(2puiq ri v ) and 



= \h)e vl ■ { Z m[H e ,P]\m)J2 A ^o 

m^h £h £m 3 

x(m\S-I j S{T-K j )\£ 1 )+ V (m|[# e ,PPi) 

x £ Au (£ 2 |S ■ I^(r - R^m) } & | (30) 



4- Hyperfine interaction 

The hyperfine interaction between the electron and nu- 
clear spins is££ 



H eI (r) 



3 



Av S ■ l 6 8{v - Kj) 



(26) 



where S = fur/2 and lj are the electron and nucleus 
spins respectively, vq = is the volume of the unit cell 
with do representing the crystal lattice parameter, and r 
(Rj) denotes the position of the electron (the j-th nu- 
cleus). A = 4/i /is/ij /(3Iv ) is the hyperfine coupling 
constant with p,Q, [Ib and fii representing the perme- 
ability of vacuum, the Bohr magneton and the nuclear 
magneton separately. 

As the Zeeman splitting of the electron is much larger 
(three orders of magnitude larger) than that of the nu- 
cleus spin, to conserve the energy for the spin relax- 
ation processes, there must be phonon-assisted transi- 
tions when considering the spin-flip processes. Tak- 
ing into account directly the BP induced motion of nu- 
clei spin of the lattice leads to a new spin relaxation 
mechanism^ 



^AttoS-Ij-MR^-VrWr- 



■R,.), (27) 



where u(R°) = VW(2^~^)(a qt) + a ql) ) e qr ) e lq ' R! ' 
is the lattice displacement vector. Therefore using the 
notation of Eq. (TTU]), $ = ^h/(2pVuj ct ^) and A ql) = 
J2j AuoS-IjV r ^(r— Rj). The second-order process of the 
surface phonon and the BP together with the hyperfine 
interaction also leads to spin relaxation: 



V 



(l 2 |<5Vc(r)|m><m|tf e7 (r)|4 



(£ 2 \H eI (r)\m)(m\5V c (r)\l 



E 



se 2 



-}<4|,(28) 



for V^P h . Similarly $ qjJ = M q7) and 



E 



(£ 2 |e iq - r |m) 



£e, — £ 
m^li 1 ••- 3 



Avo(m\S ■ lj 



S(r-R J )\£ 1 )+ V (m|e^ r Ki) 

^^o(£ 2 |S-I j( 5(r-R,)H}(4| (31) 



(3) 

for V^j_ ph . Again as the contribution from the surface 

i 2 \ 

phonon is much smaller than that of the BP, V eJ _ h can 
be neglected. It is noted that, the direct spin- phonon 
coupling due to the phonon-induced strain together with 
the hyperfine interaction gives a fourth-order scattering 
and hence induces a spin relaxation/dephasing. The in- 
teraction is 

{£ 2 \H* 8train \m){m\H eI {r)\h) 



v e f- Ph = *2>{E V " 

{l 2 \H eI {r)\m){m\H* train \l x) 



E 



Q 2 — e r. 



}(h\, (32) 




with Hg train = h z s a z /2 only changing the electron energy 
but conserving the spin polarization. It can be written 
as 



{£, yv q y - £z V q z )q z e* cl ' r ■ (33) 



Comparing this to the electron-BP interaction Eq. (|21[) . 
the ratio is « HDq/E, which is about 10 -3 . Therefore, 
the second-order term of the direct spin-phonon coupling 
due to the phonon-induced strain together with the hy- 
perfine interaction is very small and can be neglected. 
Also the coaction of the g-factor fluctuation and the hy- 
perfine interaction is very small compared to that of the 
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electron-BP interaction jointly with the hyperfine inter- 
action as hbB/E is around 10 -5 when B = 1 T. There- 
fore it can also be neglected. In the following, we only 
retain the first and the third order terms V^P_ ph and 

(3") 

V^j_ ph in calculating the spin relaxation time. 

The spin dephasing time induced by the hyperfine in- 
teraction can be calculated from the non-Markovian ki- 
netic Eq. (|18|) . for unpolarized initial nuclear spin state 
|no), resulting in 



k 



dr|^(r)| 4 cos(^|^(r)| 2 t) 



(34) 

where fk is the thermo-equilibrium distribution of the 
orbital degree of freedom. When only the lowest two 
Zeeman sublevels are considered, assuming a simple 
form of the wavefunction, |^(r)| 2 = - ^ exp(— r 2 /^) 

with du/a z representing the QD diameter/quantum well 
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width, and r 
out: 



y , the integration can be carried 



(S ( + no) ) t oc 



cos(t/to) 



sin(t/t ) 
t/t 



(35) 



Here, to = (2na z dn) / (Avo) determines the spin dephas- 
ing time. Note that to is proportional to the factor a z d 2 
where a z /d 2 is the characteristic length/area of the QD 
along z direction / in the quantum well plane. By solving 
Eq. (|18[) for various n, and summing over n, we obtain 

\\(S+)t\\ = £„ {{S+^tl- We then denne the time when 
the envelop of ||(5+)t|| decays to 1/e of its initial value 
as the spin dephasing time T<i- As mentioned above the 
hyperfine interaction can not transfer an energy of the 
order of the Zeeman splitting, thus the hyperfine inter- 
action alone can not lead to any spin relaxation; 43 

In the above discussion, the nuclear spin dipolc-dipolc 
interaction is neglected. Recently, more careful exami- 
nations based on quantum cluster expansion method or 
pair correlation method have been performedj 41 i 42 ' 43 i 47 In 
these works, the nuclear spin dipole-dipole interaction is 
also included. This interaction together with the hyper- 
fine mediated nuclear spin-spin interaction is the origin 
of the fluctuation of the nuclear spin bath. To the lowest 
order, the fluctuation is dominated by nuclear spin pair 
flipsj 41 i 42 ' 43 i 47 This fluctuation provides the source of the 
electron spin dephasing, as the electron spin is coupled 
to the nuclear spin system via hyperfine interaction. Our 
method used here includes only the hyperfine interaction 
to the second order in scattering. However, it is found 
that the dipole-dipole-interaction-induced spin dephas- 
ing is much weaker than the hyperfine interaction for a 
QD with a = 2.8 nm and do = 27 nm until the par- 
allel magnetic field is larger than ~ 20 T. 42 Therefore, 
for the situation in this paper, the nuclear dipole-dipole- 
interaction-induced spin dephasing can be ignored^ 7 - 



III. SPIN DECOHERENCE DUE TO VARIOUS 
MECHANISMS 

Following the equation-of-motion approach developed 
in Sec. II, we perform a numerical calculation of the spin 
relaxation and dephasing times in GaAs QDs. Two mag- 
netic field configurations are considered: i.e., the mag- 
netic fields perpendicular and parallel to the well plane 
(along x-axis). The temperature is taken to be T = 4 
K unless otherwise specified. For all the cases we con- 
sidered in this manuscript, the orbital level splitting is 
larger than an energy corresponding to 40 K. Therefore, 
the lowest Zeeman sublevels are mainly responsible for 
the spin decoherence. When calculating T\, the initial 
distribution is taken to be in the spin majority down 
state of the eigen-state of the Hamiltonian H e with a 
Maxwell-Boltzmann distribution /fe = Cexp[— e/ c /(fcsT)] 
for different orbital levels (C is the normalization con- 
stant). For the calculation of T2, we assign the same 
distribution between different orbital levels, but with a 
superposition of the two spin states within the same or- 
bital level. The parameters used in the calculation are 
listed in Table 

TABLE I: Parameters used in the calculation 



p 


5.3 x 10 3 kg/m 3 


K 


12.9 


Vst 


2.48 x 10 3 m/s 


9 


-0.44 


Vsl 


5.29 x 10 3 m/s 




7.0 eV 


ei4 


1.41 x 10 9 V/m 


m* 


0.067m 


A 


90 ^eV 


A33 


19.6 


70 


27.5 A 3 -eV 


I 


3 
2 


D 


1.59 x 10 4 m/s 


a 


5.6534 A 



A. Spin Relaxation Time Ti 

We now study the spin relaxation time and show how 
it changes with the well width a, the magnetic field B 
and the effective diameter do = \J hn/m*uJo. We also 
compare the relative contributions from each relaxation 
mechanism. 



1. Well width dependence 

In Fig. 1(a) and (b), the spin relaxation times induced 
by different mechanisms are plotted as function of the 
width of the quantum well in which the QD is confined 
for perpendicular magnetic field B± = 0.5 T and parallel 
magnetic field B\\ = 0.5 T respectively. We first concen- 
trate on the perpendicular magnetic field case. In Fig. 
1(a), the calculation indicates that the spin relaxation 
due to each mechanism decreases with the increase of well 
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FIG. 1: (Color online) induced by different mechanisms 
vs. the well width for (a): perpendicular magnetic field B± — 
0.5 T with (solid curves) and without (dashed curves) the 
SOC; (b) parallel magnetic field By = 0.5 T with the SOC. 
The effective diameter do = 20 nm, and temperature T — 4 K. 
Curves with ■ — Tj -1 induced by the electron-BP scattering 
together with the SOC; Curves with • — T^ 1 induced by 
the second-order process of the hyperfine interaction together 
with the BP (V e ( 7 3 l ph ); Curves with A — T' 1 induced by the 
first-order process of the hyperfine interaction together with 
the BP {V^p h ); Curves with T — T-f 1 induced by the direct 
spin-phonon coupling due to phonon-induced strain; Curves 
with ♦ — T^ 1 induced by the g- factor fluctuation. 



width. Particularly the electron-BP scattering mecha- 
nism decreases much faster than the other mechanisms. 
It is indicated in the figures that when the well width is 
small (smaller than 7 nm in the present case) , the spin re- 
laxation time is determined by the electron-BP scattering 
together with the SOC. However, for wider well widths, 
the direct spin-phonon coupling due to phonon-induced 
strain and the first-order process of hyperfine interac- 
tion combined with the electron-BP scattering becomes 
more important. The decrease of spin relaxation due to 
each mechanism is mainly caused by the decrease of the 
SOC which is proportional to a~ 2 . The SOC has two 
effects which are crucial. First, in the second order per- 
turbation the SOC contributes a finite correction to the 
Zeeman splitting which determines the absorbed/emitted 



phonon frequency and wave-vector^ Second, it leads to 
spin mixing. The decrease of the SOC thus leads to the 
decrease of Zeeman splitting and spin mixing. The for- 
mer leads to small phonon wave- vector and small phonon 
absorption/emission efficiency^ Therefore the electron- 
BP mechanism decreases rapidly with increasing a. On 
the other hand, the other two largest mechanisms can flip 
spin without the help of the SOC. The spin relaxations 
due to these two mechanisms decrease in a relatively mild 
way. It is further confirmed that without SOC they de- 
creases in a much milder way with increasing a (dashed 
curves in Fig. I). It is also noted that the spin relaxation 
rate due to the g-factor fluctuation is at least six orders 
of magnitude smaller than that due to the leading spin 
decoherence mechanisms and can therefore be neglected. 

It is noted that in the calculation, the SOC is always 
included as it has large effect on the eigen-energy and 
eigen-wavefunction of the electron^ We also show the 
spin relaxation times induced by the hyperfine interac- 
tions (V^P_ ph and Vgf_ ph ) and the direct spin-phonon 
coupling due to the phonon-induced strain but without 
the SOC as in the literature It can be seen clearly 
that the spin relaxation that includes the SOC is much 
larger than that without the SOC. For example, the spin 
relaxation induced by the second-order process of the hy- 
perfine interaction together with the BP (V^_ h ) is at 
least one order of magnitude larger when the SOC is in- 
cluded than that when the SOC is neglected. This is 
because when the SOC is neglected, {m\H e i(r)\l\) and 
(l2\H e j (r)\m) in Eq. ([29"]) are small as the matrix el- 
ements of H e i(r) between different orbital energy lev- 
els are very small. However, when the SOC is taken 
into account, the spin-up and -down levels with differ- 
ent orbital quantum numbers are mixed and therefore 
\t) and |m) include the components with the same or- 
bital quantum number. Consequently the matrix ele- 
ments of (m\H e j (r)\£i) and (t2\H e i(r)\m) become much 
larger. Therefore, spin relaxation induced by this mech- 
anism depends crucially on the SOC. 

It is emphasized from the above discussion that the 
SOC should be included in each spin relaxation mecha- 
nism. In the following calculations it is always included 
unless otherwise specified. In particular in reference to 
the mechanism of electron-BP interaction, we always con- 
sider it together with the SOC. 

We further discuss the parallel magnetic field case. In 
Fig. 1(b) the spin relaxation times due to different mech- 
anisms are plotted as function of the quantum well width 
for same parameters as Fig. 1(a), but with a parallel mag- 
netic field Bn = 0.5 T. It is noted that the spin relaxation 
rate due to each mechanism becomes much smaller for 
small a compared with the perpendicular case. Another 
feature is that the spin relaxation due to each mecha- 
nism decrease in a much slower rate with increasing a. 
The electron-BP mechanism is dominant even at a = 10 
nm but decrease faster than other mechanisms with a. It 

(3) 

is expected to be less effective than the V^ T _ h mecha- 
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nism or V^j_ ph mechanism or the direct spin-phonon cou- 
pling due to phonon-induced strain mechanism for large 
enough a. The <?-factor fluctuation mechanism is negli- 
gible again. These features can be explained as follows. 
For parallel magnetic field the contribution of the SOC to 
Zeeman splitting is much less than in the perpendicular 
magnetic field geometry.^- Moreover, this contribution is 
negative which makes Zeeman splitting smaller <2I There- 
fore, the phonon absorption/emission efficiency becomes 
much smaller for small a, i.e., large SOC. When a in- 
creases, the Zeeman splitting increases. However, the 
spin mixing decreases. The former effect is weak, and 
only cancels part of the latter, thus the spin relaxation 
due to each mechanism decrease slowly with a. 
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FIG. 3: (Color online) T^ 1 induced by different mechanisms 
vs. the parallel magnetic field B« for do = 20 nm and (a) 
a = 5 nm and (b) 10 nm. T = 4 K. Curves with ■ - 
T^ 1 induced by the electron-BP scattering; Curves with • 

— T^ 1 induced by the second-order process of the hyperfine 
interaction together with the BP (V^f_ ph ); Curves with A 

— T-^ 1 induced by the first-order process of the hyperfine 
interaction together with the BP {V^_ ph ); Curves with T 

— Tj -1 induced by the direct spin-phonon coupling due to 
phonon-induced strain; Curves with ♦ — Tj -1 induced by the 
(/-factor fluctuation. 



2. Magnetic Field Dependence 



FIG. 2: (Color online) T^ 1 induced by different mechanisms 
vs. the perpendicular magnetic field B± for do = 20 nm and 
(a) a = 5 nm and (b) 10 nm. T — 4 K. Curves with ■ — 
Tj -1 induced by the electron-BP scattering; Curves with • 

— Tj -1 induced by the second-order process of the hyperfine 
interaction together with the BP (V^p_ ph ); Curves with ▲ 

— T^ 1 induced by the first-order process of the hyperfine 
interaction together with the BP {V^j_ ph ); Curves with ▼ 

— T^ 1 induced by the direct spin-phonon coupling due to 
phonon-induced strain; Curves with ♦ — T^ 1 induced by the 
g-factor fluctuation. 



We first study the perpendicular-magnetic-field case. 
The magnetic field dependence of Ti for two different 
well widths are shown in Fig. [U^a) and Fig.[2fb). In the 
calculation, do = 20 nm. It can be seen that the ef- 
fect of each mechanism increases with the magnetic field. 
Particularly the electron-BP mechanism increases much 
faster than other ones and becomes dominant at high 
magnetic fields. For small well width (5 nm in Fig. [5^), 
the spin relaxation induced by the electron-BP scattering 
is dominant except at very low magnetic fields (0.1 T in 
the figure) where contributions from the first-order pro- 
cess of hyperfine interaction together with the electron- 
BP scattering and the direct spin-phonon coupling due to 
phonon-induced strain also contribute. It is interesting 
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to see that when a is increased to 10 nm, the electron-BP 
scattering is the largest spin relaxation mechanism only 
at high magnetic fields (>1.1 T). For 0.4 T < B± < 1.1 
T (Bj_ < 0.4 T), the direct spin-phonon coupling due 
to the phonon-induced strain (the first order hyperfine 
interaction together with the BP ) becomes the largest 
relaxation mechanism. It is also noted that there is no 
single mechanism which dominates the whole spin relax- 
ation. Two or three mechanisms are jointly responsible 
for the spin relaxation. It is indicated that the spin relax- 
ations induced by different mechanisms all increase with 
B±. This can be understood from a perturbation theory: 
when the magnetic field is small the spin relaxation be- 
tween two Zeeman split states for each mechanism is pro- 
portional to n(AE)(AE) m (AE is the Zeeman splitting) 
with ttl = 7 for electron-BP scattering due to the de- 
formation potential 18 "^ and for the second-order process 
of the hyperfine interaction together with the electron- 
BP scattering due to the deformation potential V^f}_ ph j22 
rn = 5 for electron-BP scattering due to the piezoelec- 
tric couplin g 15 i 18 i 25 and for the second-order process of 
the hyperfine interaction together with the electron-BP 
scattering due to the piezoelectric coupling V^_ ph -^ and 
777. = 5 for the direct spin-phonon coupling due to phonon- 
induced strain^ m = 1 for the first-order process of the 
hyperfine interaction together with the BP V^'_ ph . The 
spin relaxation induced by the ^-factor fluctuation is pro- 
portional to n(AE)(AE) 5 B\. For most of the cases stud- 
ied, AE is smaller than k B T, hence n(AE) ~ k B T/AE, 
andn(AE)(AE) m ~ (AE)" 1 - 1 . m > 1 hold for all mech- 
anism except the V^P_ ph mechanism, therefore the spin 
relaxation due to these mechanisms increases with in- 
creasing B±. However, from Eq. (f2"T)) one can see that 
it has a term with V r , which indicates that the effect of 
this mechanism is proportional to I/do- As the vector 
potential of the magnetic field increases the confinement 
of the QD and gives rise to smaller effective diameter do, 
this mechanism also increases with the magnetic field in 
the perpendicular magnetic field geometry. 

We then study the case with the magnetic field par- 
allel to the quantum well plane. In Fig. [3] the spin re- 
laxation induced by different mechanisms are plotted as 
function of the parallel magnetic field Bn for two dif- 
ferent well widths. In the calculation, do = 20 nm. It 
can be seen that, similar to the case with perpendicular 
magnetic field, the effects of most mechanisms increase 
with the magnetic field. Also the electron-BP mechanism 
increases much faster than the other ones and becomes 
dominant at high magnetic fields. However, without the 
orbital effect of the magnetic field in the present con- 
figuration, the effect of h changes very little with 
the magnetic field. For both small (5 nm in Fig. [3ja)) 
and large (10 nm in Fig. [3Jb)) well widths, the electron- 
BP scattering is dominant except at very low magnetic 
field (0.1 T in the figure), where the first-order process of 
the hyperfine interaction together with the electron-BP 
interaction V^p_ , also contributes. 




d Q (nm) 

FIG. 4: (Color online) T-j -1 induced by different mechanisms 
vs. the effective diameter do for B± — 0.5 T and (a) a = 5 nm 
and (b) 10 nm. T = 4 K. Curves with ■ — Tf 1 induced by 
the electron-BP scattering; Curves with • - T x _ induced by 
the second-order process of the hyperfine interaction together 
with the BP (V}fl ph ); Curves with A — Tf 1 induced by the 
first-order process of the hyperfine interaction together with 
the BP (Vjjl , ); Curves with T — Tf 1 induced by the direct 
spin-phonon coupling due to phonon-induced strain; Curves 
with ♦ — Tj -1 induced by the gr-factor fluctuation. 



3. Diameter Dependence 

We now turn to the investigation of the diameter de- 
pendence of the spin relaxation. We first concentrate on 
the perpendicular magnetic field geometry. The spin re- 
laxation rate due to each mechanism is shown in Fig. 0^ 
for a small (a — 5 nm) and Fig. [4Jd for a large (a = 10 
nm) well widths respectively with a fixed perpendicu- 
lar magnetic field B± = 0.5 T. In the figure, the spin 
relaxation rate due each mechanism except V^_ ph in- 
creases with the effective diameter. Specifically, the ef- 
fect of the electron-BP mechanism increases very fast, 
while the effect of the direct spin-phonon coupling due to 
phonon-induced strain mechanism increases very mildly. 
The V^jl_ h decreases with do slowly. Other mechanisms 
are unimportant. The electron-BP mechanism eventu- 
ally dominates spin relaxation when the diameter is large 
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laxation due to V^_ ph decreases with the diameter. This 

is because h contains a term V r [Eq. ([2T)l] which 

decreases with the increase of do- Physically speaking, 
the decrease of the effect of V^'_ ph is due to the fact 
that the spin mixing due to the hyperfine interaction de- 
creases with the increase of the number of nuclei within 
the dot N as the random Overhauser field is proportional 
to 1/y/W. The spin relaxation induced by the ^-factor is 
also negligible here for both small and large well width. 

We then turn to the parallel magnetic field case. In 
the calculation, Sm = 0.5 T. The results are shown for 
both small well width (a = 5 nm in Fig.[5][a)) and large 
well width (a = 10 nm in Fig. EJb)) respectively. Simi- 
lar to the perpendicular magnetic field case, the effect of 
every mechanism except the V^j_ h mechanism increases 
with increasing diameter. The effect of the electron-BP 
mechanism increases fastest and becomes dominant for 
do > 12 nm for both small and large well width. For 
do < 12 nm for the two cases the first-order process 
of the V^P_ ph mechanism becomes dominant. The ef- 

(3) 

feet of the V^_ ph mechanism become larger than that of 
the direct spin-phonon coupling due to phonon-induced 
strain mechanism. However, these two mechanism are 
still unimportant and becomes more and more unimpor- 
tant for larger do- Here, the spin relaxation induced by 
the g-factor is negligible. 



FIG. 5: (Color online) T 1 1 induced by different mechanisms 
vs. the effect diameter do with Bu = 0.5 T and (a) a = 5 nm 
and (b) 10 nm. T = 4 K. Curves with ■ — Tf 1 induced by 
the electron-BP scattering; Curves with • — T^ 1 induced by 
the second-order process of the hyperfine interaction together 
with the BP (V^l ph ); Curves with A — T!f 1 induced by the 
first-order process of the hyperfine interaction together with 
the BP (V}jl h ); Curves with T — T{~ 1 induced by the direct 
spin-phonon coupling due to phonon-induced strain; Curves 
with ♦ — Tj - induced by the g-factor fluctuation. 



enough. The threshold increases from 12 nm to 26 nm 
when the well width increases from 5 nm to 10 nm. For 
small diameter the V^j_ h and the direct spin-phonon 
coupling due to phonon-induced strain mechanism dom- 
inate the spin relaxation. The increase/decrease of the 
spin relaxation due to these mechanisms can be under- 
stood from the following. The effect of the SOC on the 
Zeeman splitting is proportional to dp f° r small magnetic 
field^ The increase of do thus leads to a increase of Zee- 
man splitting, therefore the efficiency of the phonon ab- 
sorption/emission increases. Another effect is that the in- 
crease of do will increase the phonon absorption/emission 
efficiency due to the increase of the form factor Thus 
the spin relaxation increases. Moreover, the spin mix- 
ing is also proportional to do alsoJ^ This leads to much 
faster increasing of the effect of the electron-BP mecha- 

(3) 

nism and the h mechanism. However, the spin re- 



4- Comparison with Experiment 

In this subsection, we apply our analysis to experiment 
data in Ref. 0. We first show that our calculation is in 
good agreement with the experimental results. Then we 
compare contributions from different mechanisms to spin 
relaxation as function of the magnetic field. In the cal- 
culation we choose the quantum dot diameter do = 56 
nm (Hu>o = 1-1 meV as in experiment). The quantum 
well is taken to be an infinite-depth well with a = 13 
nm. The Dresselhaus SOC parameter 7o(fcz) i s taken 
to be 4.5 meV-Aand the Rashba SOC parameter is 3.3 
meV-A. T — K as ksT <C g/J-BB in the experiment. 
The magnetic field is applied parallel to the well plane in 
[110]-direction. The Dresselhaus cubic term is also taken 
into consideration. All these parameters are the same 
with (or close to) those used in Ref. [24| in which a cal- 
culation based on the electron-BP scattering mechanism 
agrees well with the experimental results. For this mech- 
anism, we reproduce their results. The spin relaxation 
time measured by the experiments (black dots with er- 
ror bar in the figure) almost coincide with the calculated 
spin relaxation time due to the electron-BP scattering 
mechanism (curves with ■ in the figure)^ It is noted 
from the figure that other mechanisms are unimportant 
for small magnetic field. However, for large magnetic 
field the effect of the direct-spin phonon coupling due 
to phonon-induced strain becomes comparable with that 
of the electron-BP mechanism. At B\\ — 10 T, the two 
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differs by a factor of ~ 5. 
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FIG. 6: (Color online) T^ 1 induced by different mechanisms 
vs. the parallel magnetic field Bu in the [110] direction for 
do = 56 nm and a = 13 nm with both the Rashba and Dres- 
selhaus SOCs. T = K. The black dots with error bar is the 
experimental results in Ref.|7|. Curves with ■ — Tj -1 induced 
by the electron-BP scattering; Curves with • — Tj -1 induced 
by the second-order process of the hyperfine interaction to- 
gether with the BP (V}p_ ph ); Curves with A — Tf 1 induced 
by the first-order process of the hyperfine interaction together 
with the BP (V}jl ph ); Curves with T — Tf 1 induced by the 
direct spin-phonon coupling due to phonon- induced strain; 
Curves with ♦ — T{~ induced by the g-factor fluctuation. 



B. Spin Dephasing Time Ti 

In this subsection, we investigate the spin dephasing 
time for different well widths, magnetic fields and QD 
diameters. As in the previous subsection, the contribu- 
tions of the different mechanisms to spin dephasing are 
compared^ To justify the first Born approximation in 
studying the hyperfine interaction induced spin dephas- 
ing, we focus mainly on the high magnetic field regime 
of B > 3.5 T. A typical magnetic field is 4 T. We also 
demonstrate via extrapolation that in the low magnetic 
field regime spin dephasing is dominated by the hyperfine 
interaction. 



1. Well Width Dependence 

In Fig.[7Jthe well width dependence of the spin dephas- 
ing induced by different mechanisms is presented under 
the perpendicular (a) and parallel (b) magnetic fields. In 
the calculations B± = 4 T/By = 4 T and d = 20 nm. It 
can be seen in both figures that the spin dephasing due to 
each mechanism decreases with a. Moreover, the spin de- 
phasing due to the electron-BP scattering decreases much 
faster than that due to the hyperfine interaction. These 
features can be understood as following. The spin de- 
phasing due to electron-BP scattering depends crucially 



on the SOC. As the SOC is proportional to a -2 , the spin 
dephasing decreases fast with a. For the hyperfine inter- 
action, from Eq. ([33]) one can deduce that the decay rate 
of ||(5+)t|| is mainly determined by the factor l/(a z d|) 
(here a z = a), which thus decreases with a, but in a very 
mild way. The fast decrease of the electron-BP mecha- 
nism makes it eventually unimportant. For the present 
perpendicular-magnetic-field case the threshold is around 
2 nm. For parallel magnetic field it is even smaller. A 
higher temperature may enhance the electron-BP mech- 
anism (see discussion in Sec. V) and make it more im- 
portant than the hyperfine mechanism. It is noted that 
other mechanisms contribute very little to the spin de- 
phasing. Thus, in the following discussion, we do not con- 
sider these mechanisms. Comparing Figs UJa) and (b), 
one finds that a main difference is that the electron-BP 
mechanism is less effective for the parallel-magnetic-field 
case. As has been discussed in the previous subsection, 
the spin mixing and the Zeeman splitting in the parallel 
filed case is smaller than those in the perpendicular field 
case. Therefore, the electron-BP mechanism is weakened 
markedly. 

Similar to Fig. [TJ the SOC is always included in the 
computation as it has large effect on the eigen-energy and 
eigen-wavefunction of the electrons. The spin dephasings 
calculated without the SOC for the hyperfine interaction, 
the direct spin-phonon coupling due to phonon-induced 
strain and the g-factor fluctuation are also shown in Fig. 
[7£a) as dashed curves. It can be seen from the figure 
that for the spin dephasings induced by the direct spin- 
phonon coupling due to phonon-induced strain and by 
the g-factor fluctuation, the contributions with the SOC 
are much larger than those without. This is because when 
the SOC is included, the fluctuation of the effective field 
induced by both mechanisms becomes much stronger and 
more scattering channels are opened. However, what 
should be emphasized is that the spin dephasings induced 
by the hyperfine interaction with and without the SOC 
are nearly the same (the solid and the dashed curves 
nearly coincide). That is because the change of the wave- 
function ^(r) due to the SOC is very small (less than 1 
% in our condition) and therefore the factor l/(a z cZ 2 ) is 
almost unchanged when the SOC is neglected. Thus the 
spin dephasing rate is almost unchanged. 

In the inset of Fig. [7(a), the time evolution of 1 1 (S + ) t \ \ 
induced by the hyperfine interaction is shown, with a = 2 
nm. It can be seen that 1 1 (S+) t 1 1 decays very fast and de- 
creases to less than 10 % of its initial value within the 
first two oscillating periods. Therefore, Ti is determined 
by the first two or three periods of ||(5+)t||. Thus the 
correction of the long time dynamics due to higher or- 
der scattering 33 contributes little to the spin dephasing 
time. For quantum computation and quantum informa- 
tion processing, the initial, e.g., 1 % decay of ||(5+)t|| 
may be more important than the 1/e deca y 42 i 43 Indeed, 
the spin dephasing time defined by the exponential fitting 
of 1 % decay is short than that defined by the 1/e decay. 
However, the two differs less than 5 times. For a rough 
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FIG. 7: (Color online) T 2 1 induced by different mechanisms 
vs. the well width for d = 20 nm. T = 4 K. (a): B± = 4 T 
with (solid curves) and without (dashed curves) the SOC; (b): 
By = 4 T only with the SOC. Curve with ■ — T 2 ~ x induced 
by the electron-BP interaction; Curves with • — T 2 _1 induced 
by the hyperfine interaction; Curves with T — T 2 _1 induced 
by the direct spin-phonon coupling due to phonon-induced 
strain; Curves with ♦ — T 2 ~ induced by g-factor fluctuation; 
A — T 2 _1 induced by the second-order process of the hyper- 
fine interaction together with the BP (V}f}_ ph ); Curves with □ 
— T 2 _1 induced by the first-order process of the hyperfine in- 
teraction together with the BP {V^_ ph ). The time evolution 
of ||(iS+)t|| induced by the hyperfine interaction with a = 2 
nm is shown in the inset of (a). 



comparison of contributions from different mechanisms to 
spin dcphasing where only the order-of-magnitude differ- 
ence is concerned (see Figs. 7-9), this difference due to 
the definition does not jeopardize our conclusions. 



2. Magnetic Field Dependence 

We then investigate the magnetic field dependence of 
the spin dephasing induced by the electron-BP scatter- 
ing and by the hyperfine interaction for two different well 
widths (a = 3 nm and a = 5 nm) with both perpendicular 
and parallel magnetic field. From Fig. |3(a) and (b) one 



FIG. 8: (Color online) T 2 ~ x induced by the electron-BP scat- 
tering and the hyperfine interaction vs. (a): the perpendic- 
ular magnetic field B± ; (b): the parallel magnetic field Bu 
for a = 3 nm (solid curves) and 5 nm (dashed curves). T — 4 
K, and do = 20 nm. Curves with ■ — T 2 induced by the 
electron-BP interaction; Curves with • — T 2 induced by the 
hyperfine interaction. 



can see that the spin dephasing due to the electron-BP 
scattering increases with magnetic field, whereas that due 
to the hyperfine interaction decreases with magnetic field. 
Thus, the electron-BP mechanism eventually dominates 
the spin dephasing for high enough magnetic field. The 
threshold is B° ± = 4 T / B?, = 7 T for a = 3 nm with per- 
pendicular/parallel magnetic field. For larger well width, 
e.g., a — 5 nm with parallel magnetic field or perpendicu- 
lar magnetic field, the threshold magnetic fields increase 
to larger than 8 T. The different magnetic field depen- 
dences above can be understood as following. Besides 
spin relaxation, the spin-flip scattering also contributes 
to spin dephasingi^i As has been demonstrated in Sec. 
1IIA, the electron-BP scattering induced spin-flip tran- 
sition rate increases with the magnetic field. Therefore 
the spin dephasing rate increases with the magnetic field 
also. In contrast, spin dephasing induced by the hyper- 
fine interaction decreases with the magnetic field. This 
is because when the magnetic field becomes larger, the 
fluctuation of the effective magnetic field due to the sur- 
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rounding nuclei becomes insignificant compared. There- 
fore, the hyperfine-interaction-induced spin dephasing is 
reduced. Similar results have been obtained by Deng and 
Hu4i 
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FIG. 9: (Color online) T 2 1 induced by the electron-BP scat- 
tering and the hyperfine interaction vs. the effective diameter 
do T = 4 K. (a): B ± = 4 T ; (b): By = 4 T for a = 3 nm 
(solid curves) and 5 nm (dashed curves). Curves with ■ — 
T 2 -1 induced by the electron-BP interaction; Curves with • 
— T 2 _1 induced by the hyperfine interaction. 



3. Diameter Dependence 

In Fig. [5] the spin dephasing times induced by the 
electron-BP scattering and the hyperfine interaction are 
plotted as function of the diameter do for a small (a = 3 
nm) and a large (a = 5 nm) well widths. In the cal- 
culation, B x = 4 T in Fig. [^a) and By = 4 T in (b). 
It is noted that the effect of the electron-BP mechanism 
increases rapidly with do, whereas the effect of the hy- 
perfine mechanism decreases slowly. Consequently, the 
electron-BP mechanism eventually dominates the spin 
dephasing for large enough do- The threshold is dg = 19 
(27) nm for a = 3 (5) nm case with the perpendicular 
magnetic field and d§ = 26 (30) nm for a = 3 (5) nm 
case under the parallel magnetic field. As has been dis- 



cussed in Sec. IIIA, both the effect of the SOC and the ef- 
ficiency of the phonon absorption/emission increase with 
do- Therefore, the spin dephasing due to the electron-BP 
mechanism increases rapidly with do' 18 ' 21 The decrease 
of the effect of the hyperfine interaction is due to the de- 
crease of the factor l/(a z dT) [Eq. [33] with the diameter 
d . 



IV. SPIN RELAXATION TIMES FROM FERMI 
GOLDEN RULE AND FROM EQUATION OF 
MOTION 

In this section, we will try to find a proper method to 
average over the transition rates from the Fermi Golden 
rule, Tjl+f; to give the spin relaxation time T%. In the 
limit of small SOC, we rederive Eq. ([T]) from the equation 
of motion. We further show that Eq. (Q} fails for large 
SOC where a full calculation from the equation of motion 
is needed. 
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FIG. 10: (Color online) Spin relaxation time Ti calculated 
from the equation-of-motion approach (■) v.s. that obtained 
from Eq. (1) (•) as function of (a): the strength of the SOC 
for T=12 K; (b): the temperature for 7 = 70. The well width 
a = 5 nm, perpendicular magnetic field B± = 0.5 T, QD 
diameter do = 30 nm. The ratio of the two R is also plotted 
in the figure. Note the scale of T^ 1 is at the right hand side 
of the frame. 
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We first rederive Eq. ([T]) for small SOC from the equa- 
tion of motion. In QDs, the orbital level splitting is 
usually much larger than the Zeeman splitting. Each 
Zeeman sublevel has two states: one with majority up- 
spin, the other with majority down-spin. We call the 
former "minus state" (as it corresponds to a lower en- 
ergy) while the latter "plus state". For small SOC, 
the spin mixing is small. Thus we neglect the much 
smaller contribution from the off-diagonal terms of the 
density matrix to S z . Therefore S z (t) = fi±(t) 
where i± denotes the plus/minus state of the z-th or- 
bital state. For small SOC, the spin relaxation is 
much slower than the orbital relaxation i 25 i 55 This im- 
plies that the time takes to establish equilibrium within 
the plus/minus states is much smaller than the spin relax- 
ation time. Thus we can assume a equilibrium (Maxwell- 
Boltzmann) distribution between the plus/minus states 
at any time. The distribution function is therefore given 
by fi±{t) = N ± (t)exp(-e l± /k B T)/Z ± . Here N ± (t) = 
J2i /i±(*) i s the total probability of the plus/minus states 
with N + (t) + N_(t) = 1 for single electron in QD and 
Z± = '^2 i Gxp(—ei±/kBT) is the partition function for 
the plus/minus state. At equilibrium, N± = N± . The 
equation for S z (t) is hence, 

!*.(*) = j t [sS)-s?] 

= J2 S l ±e M-^±/k B T)/Z ± j t SN ± (t) 7 (36) 

i± 

with 5N±(t) = N±(t) ~ N"^. As the orbital level 
splitting is usually much larger than the Zeeman split- 
ting, the factor exp(—Si±/kBT)/Z± can be approximated 
by ex_p(-e io /k B T)/Zo with e i0 = \{e l+ + £j_) and 
Zq = J2i ex P[ — Sio/kBT]. Further using the particle- 
conservation relation J2± 5N±(t) = 0, one has 

j f S z (t) = [J2(Si + -Si-)eM-^o/kBT)/Z ]j t 6N + (t) . 

(37) 

As S z (t) - S? = [8N + {t)/Z ]Ei(Sp 
SI ) exp(— Sio/hsT), one finds that the spin relax- 
ation time is nothing but the relaxation time of JV+. 
The next step is to derive the equation of 4:SN+(t), 
which is given in our previous work 4^ 

i 

= -Yi^f-+r^ f+ ] e ~ ei °^ BT 6N + (t) .(38) 
»>/ 

Thus spin relaxation time is given by, 

1 «^ p -e,a/k B T 

¥l = T,^f- + 7 T- 1 -. /+ )^ r - • 09) 



Furthermore, substituting e £i °/ kBT /Zo by ff± = 
exp(— Ei±/kBT)/Z±, we have 

Y = £^+-/-/°+ + r-}^ f+ fl) ■ (40) 
1 i,f 

This is exactly Eq. (|TJ) . 
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FIG. 11: (Color online) Spin relaxation time Ti, spin dephas- 
ing time T2 and T1/T2 against temperature T. B± = 4 T, 
a = 5 nm and do = 30 nm. Note the scale of Ti and T2 is at 
the right hand side of the frame. 

For large SOC, or large spin mixing due to anticrossing 
of different spin states ; 19 i 25 the spin relaxation rate be- 
comes comparable with the orbital relaxation rate. Fur- 
thermore, the decay of the off-diagonal term of the den- 
sity matrix should contribute to the decay of S z . There- 
fore, the above analysis does not hold. In this case, it is 
difficult to obtain such a formula, and a full calculation 
from the equation-of-motion is needed. 

In Fig. HUIa), we show that (for T = 12 K, a = 5 
nm, B± = 0.5 T, do = 30 nm) the spin relaxation times 
T\ calculated from equation-of-motion approach and that 
obtained from Eq. (|40[) . Here, for simplicity and without 
loss of generality, we consider only the electron-BP scat- 
tering mechanism. The discrepancy of T\ obtained from 
the two approaches increases with 7. At 7 = IO70, the 
ratio of the two becomes as large as ~ 3. In Fig. fTUTb), 
we plot the spin relaxation times obtained via the two 
approaches as function of temperature for 7 = 70 with 
other parameters remaining unchanged. It is noted that 
the discrepancy of T± obtained from the two approaches 
increases with temperature. For high temperature, the 
higher levels are involved in the spin dynamics where the 
SOC becomes larger. At 40 K, the discrepancy is as large 
as 60 %. The ratio increases very slowly for T < 20 K 
where only the lowest two Zeeman sublevels are involved 
in the dynamics. 
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FIG. 12: (Color online) Spin relaxation time Ti against tem- 
perature T for (a): B± = 0.5 T; (b): B± = 0.9 T. a = 10 
nm and do = 20 nm. Curves with ■ — T[~ induced by the 
electron-BP scattering together with the SOC; Curves with • 
— T^ 1 induced by the second-order process of the hyperfine 
interaction together with the BP (V^_ ph )\ Curves with A — 
T^ 1 induced by the first-order process of the hyperfine inter- 
action together with the BP (V^j_ ph ); Curves with T — Tj -1 
induced by the direct spin-phonon coupling due to phonon- 
induced strain; Curves with ♦ — 77/ 1 induced by the g- factor 
fluctuation. 



V. TEMPERATURE DEPENDENCE OF SPIN 
RELAXATION TIME Ti AND SPIN DEPHASING 
TIME T 2 

We first study the relative magnitude of the spin re- 
laxation time T\ and the spin dephasing time T2. We 
consider a QD with do = 30 nm and a = 5 nm at B± = 4 
T where the largest contribution to both spin relaxation 
and dephasing comes from the electron-BP scattering 
(see Fig. 0|a) and Fig. [5{a), we have checked that the 
electron-BP scattering mechanism is dominant through- 
out the temperature range). From Fig.QTJ one finds that 
when the temperature is low (T < 5 K in the figure), 
Ti = 2T±, which is in agreement with the discussion in 
Ref. However, T1/T2 increases very quickly with T 
and for T = 20 K, Ti/T 2 - 2 x 10 2 . This is understood 



from the fact that when T is low, the electron mostly dis- 
tributes in the lowest two Zeeman sublevels. For small 
SOC, Golovach et al. have shown via perturbation the- 
ory that phonon induces only the spin-flip noise in the 
leading order. Consequently, T2 = 2T\*& When the tem- 
perature becomes comparable with the orbital level split- 
ting Hluq, the distribution over the upper orbital levels is 
not negligible any more. As mentioned previously, the 
SOC contributes a non-trivial part to the Zeeman split- 
ting. Specifically, the second order energy correction due 
to the SOC contributes to the Zeeman splitting. The 
energy correction for different orbital levels is generally 
unequal (always larger for higher levels). When the elec- 
tron is scattered by phonons randomly from one orbital 
state to another one with the same major spin polariza- 
tion, the frequency of its precession around z direction 
changes. Continuous scattering leads to random fluctu- 
ation of the precession frequency and thus leads to spin 
dephasing,— ^ Note that this fluctuation only leads to a 
phase randomization of S+, but not flips the z compo- 
nent spin S z , i.e., not leads to spin relaxation. There- 
fore, the spin dephasing becomes stronger than the spin 
relaxation for high temperatures. Moreover, this effect 
increases with temperature rapidly as the distribution 
over higher levels and the phonon numbers both increase 
with temperature. 

We further study the temperature dependence of spin 
relaxation for lower magnetic field and larger quantum 
well width where other mechanisms may be more im- 
portant than the electron-BP mechanism. In Fig. [T27 a). 
the spin relaxation time is plotted as function of tem- 
perature for B± — 0.5 T, a = 10 nm and c?o — 20 nm. 
It is seen from the figure that the direct spin-phonon 
coupling due to phonon-induced strain mechanism dom- 
inates the spin relaxation throughout the temperature 
range. It is also noted that for T < 4 K the spin relax- 
ation rates induced by different mechanisms all increase 
with temperature according to the phonon number fac- 
tor 2n(E z i) + 1 with E z \ being the Zeeman splitting of 
the lowest Zeeman sublevels. However, for T > 4 K, the 
spin relaxation rates induced by the direct spin-phonon 
coupling due to phonon-induced strain and the electron- 
BP interaction increase rapidly with temperature, while 
the spin relaxation rates induced by V^_ ph and V^P h 
increase mildly according to 2n{E z \) + 1 throughout the 
temperature range. These features can be understood as 
what follows. For T < 4 K, the distribution over the high 
levels is negligible. Only the lowest two Zeeman sublevels 
involve in the spin dynamics. The spin relaxation rates 
thus increase with 2n[E z \) + 1 and the relative impor- 
tance of each mechanism does not change. Therefore, our 
previous analysis on comparison of relative importance of 
different spin decoherence mechanisms at 4 K holds true 
for the range < T < 4 K. When the temperature gets 
higher, the contribution from higher levels becomes more 
important. Although the distribution at the higher lev- 
els is still very small, for the direct spin-phonon coupling 
mechanism, the transition rates between the higher lev- 
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els and that between higher levels and the lowest two 
sublevels are very large. For the electron-BP mechanism 
the transition rates between the higher levels are very 
large due to the large SOC in these levels. Therefore, the 
contribution from the higher levels becomes larger than 
that from the lowest two sublevels. Consequently, the in- 
crease of temperature leads to rapid increase of the spin 
relaxation rates. However, for the two hyperfine mech- 
anisms: the V^j_ ph and the V^f!_ ph , the spin relaxation 
rates does not change much when the higher levels are 
involved. They thus increase by the phonon number fac- 
tor. 

In Fig. [T27 b) we show the temperature dependence of 
the spin relaxation time for the same condition but with 
B = 0.9 T. It is noted that the spin relaxation rate due 
to the electron-BP mechanism catches up with that in- 
duced by the direct spin-phonon coupling due to phonon- 
induced strain at T — 9 K and becomes larger for higher 
temperature. This indicates that the temperature depen- 
dence of the two mechanisms are quite different. 

In Fig. [13] we show the spin dephasing induced by 
electron-BP scattering and the hyperfine interaction as 
function of temperature for B± = 4 T, a = 10 nm and 
do = 20 nm. We choose the conditions so that the spin 
dephasing is dominated by the hyperfine interaction at 
low temperature. However, the effect of the electron-BP 
mechanism increases with temperature quickly while that 
of the hyperfine interaction remains nearly unchanged. 
The fast increase of the effect from the electron-BP scat- 
tering is due to three factors: 1) the increase of the 
phonon number; 2) the increase of scattering channels; 
and 3) the increase of the SOC induced spin mixing in 
higher levels. On the other hand, from Eq. [351 one can 
deduce that the spin dephasing rate of the hyperfine in- 
teraction depends mainly on the factor l/(a z d?,) with 

a z /d?< is the characteristic length/area along the z di- 
rection / in the quantum well plane. For higher levels, 
the d?< is larger, but only about a factor smaller than 10. 
Thus the effect of the hyperfine interaction increases very 
slowly with temperature. 

It should be noted that in the above discussion, we 
neglected the two-phonon scattering mechanism ) 15 ' 46 i 50 
which may be important at high temperature. The con- 
tribution of this mechanism should be calculated via 
the equation-of-motion approach developed in this paper, 
and compared with the contribution of other mechanisms 
showed here. 



VI. CONCLUSION 

In conclusion, we have investigated the longitudinal 
and transversal spin decoherence times Ti and T 2 , called 
spin relaxation time and spin dephasing time, in differ- 
ent conditions in GaAs QDs from the equation-of-motion 
approach. Various mechanisms, including the electron- 
BP scattering, the hyperfine interaction, the direct spin- 
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FIG. 13: (Color online) Spin relaxation time Ti against tem- 
perature T. B± = 4 T, a — 10 nm and do = 20 nm. Curves 
with ■ — T 2 _1 induced by the electron-BP scattering together 
with the SOC; Curves with • — T 2 _1 induced by the hyperfine 
interaction 



phonon coupling due to phonon-induced strain and the 
g-factor fluctuation are considered. Their relative im- 
portance is compared. There is no doubt that for spin 
decoherence induced by electron-BP scattering, the SOC 
must be included. However, for spin decoherence induced 
by the hyperfine interaction, the direct spin-phonon cou- 
pling due to phonon-induced strain, g-factor fluctua- 
tion, and hyperfine interaction combined with electron- 
phonon scattering, the SOC is neglected in the existing 
literaturei 27 ' 28 ! 45 Our calculations have shown that, as 
the SOC has marked effect on the eigen-energy and the 
eigen-wavefunction of the electron, the spin decoherence 
induced by these mechanisms with the SOC is larger than 
that without it. Especially, the decoherence from the 
second-order process of hyperfine interaction combined 
with the electron-BP interaction increases at least one 
order of magnitude when the SOC is included. Our cal- 
culations show that, with the SOC, in some conditions 
some of these mechanisms (except (/-factor fluctuation 
mechanism) can even dominate the spin decoherence. 

There is no single mechanism which dominates spin re- 
laxation or spin dephasing in all parameter regimes. The 
relative importance of each mechanism varies with the 
well width, magnetic field and QD diameter. In particu- 
lar, the electron-BP scattering mechanism has the largest 
contribution to spin relaxation and spin dephasing for 
small well width and/or high magnetic field and/or large 
QD diameter. However, for other parameters the hyper- 
fine interaction, the first-order process of the hyperfine 
interaction combined with electron-BP scattering, and 
the direct spin-phonon coupling due to phonon-induced 
strain can be more important. It is noted that the <?- 
factor fluctuation always has very little contribution to 
spin relaxation and spin dephasing which can thus be 
neglected all the time. For spin dephasing, the electron- 
BP scattering mechanism and the hyperfine interaction 
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mechanism are more important than other mechanisms 
for magnetic field higher than 3.5 T. For this regime, 
other mechanisms can thus be neglected. It is also shown 
that spin dephasing induced by the electron-BP mecha- 
nism increases rapidly with temperature. Extrapolated 
from our calculation, the hyperfine interaction mecha- 
nism is believed to be dominant for small magnetic held. 

We also discussed the problem of finding a proper 
method to average over the transition rates t~^j obtained 
from the Fermi Golden rule, to give the spin relaxation 
time T\ at finite temperature. For small SOC, we red- 
erived the formula for T\ at finite temperature used in 
the existing literatur o 18 i 51 ' 52 from the equation of mo- 
tion. We further demonstrated that this formula is in- 
adequate at high temperature and/or for large SOC. For 
such full calculation from the equation-of-motion 

approach is needed. The equation-of-motion approach 
provides an easy and powerful way to calculate the spin 
decoherence at any temperature and SOC. 



We also studied the temperature dependence of spin re- 
laxation Ti and dephasing T 2 - We show that for very low 
temperature if the electron only distributes on the low- 
est two Zeeman sublevels, T 2 = 2T\. However, for higher 
temperatures, the electron spin dephasing increases with 
temperature much faster than the spin relaxation. Con- 
sequently T\ 3> T%. The spin relaxation and dephasing 
due to different mechanisms are also compared. 
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